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Abstract 

A measurement is presented of the charged hadron multiplicity in hadronic PbPb 
collisions, as a function of pseudorapidity and centrality, at a collision energy of 
2.76 TeV per nucleon pair. The data sample is collected using the CMS detector and a 
minimum-bias trigger, with the CMS solenoid off. The number of charged hadrons is 
measured both by counting the number of reconstructed particle hits and by forming 
hit doublets of pairs of layers in the pixel detector. The two methods give consistent 
results. The charged hadron multiplicity density, dNch/drj\rj^o, for head-on collisions 
is found to be 1612 ± 55, where the uncertainty is dominated by systematic effects. 
Comparisons of these results to previous measurements and to various models are 
also presented. 
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1 Introduction 

Quantum chromodynamics (QCD), the theory of strong interactions, predicts a phase transi- 
tion at high temperature between hadronic and deconfined matter 111. Strongly interacting 
matter under extreme conditions can be studied experimentally using ultrarelativistic colli- 
sions of heavy nuclei. The field entered a new era in November 2010 when the Large Hadron 
Collider (LHC) produced the first PbPb collisions at a centre-of-mass energy per nucleon pair 
of 2.76 TeV. This represents an increase of more than one order of magnitude over the highest- 
energy nuclear collisions previously achieved in the laboratory. The multiplicity of charged 
particles produced in the central-rapidity region is a key observable characterising the proper- 
ties of the quark-gluon matter created in these collisions [2J. 

Nuclei are extended objects, and their collisions occur at various impact parameters, referred 
to as "centralities". The studies of the dependence of the charged particle density on the type 
of colliding nuclei, on the centre-of-mass energy, and on the collision geometry are important 
for understanding the relative contributions of hard scattering and soft processes to particle 
production and provide insight into the partonic structure of the nuclei. 

In this paper we report measurements of the multiplicity density dN^h/d?/ of primary charged 
hadrons. The analysis is based on the 2.76 TeV-per-nucleon PbPb collision data recorded by 
the Compact Muon Solenoid (CMS) detector in December 2010, in runs without magnetic field. 
The pseudorapidity is defined as // = — ln[tan(0/2)] with 9 the polar angle with respect to the 
counterclockwise beam direction (the z axis). The number of primary charged hadrons N^h is 
defined as all charged hadrons produced in an event including decay products of particles with 
proper lifetimes less than 1 cm. 

A detailed description of the CMS experiment can be found in Ref. HSj. The pixel tracker used 
for the analysis covers the region \tj\ < 2.5 and a full 2nin azimuth, with 66M detector channels 
out of which 97.5% were functional during data taking. It consists of a three-layer barrel pixel 
detector (BPIX) and two endcap disks at each barrel end. Only the barrel section was used 
in this analysis. The first BPIX layer is located at a radius betw^een 3.6 and 5.2 cm from the 
beam line, the second between 6.6 and 8.0 cm, and the third between 9.4 and 10.8 cm. The 
detectors used for event selection are the hadron forward (HP) calorimeters, w^hich cover the 
region 2.9 < |;/| < 5.2, the beam scintillator counters (BSC), in the range 3.23 < \ri\ < 4.65, 
and the beam pick-up timing (BPTX) devices located at z = ±176 m from the interaction point. 
The operation of the CMS detector w^ith zero magnetic field has the benefit of an increase in 
the acceptance for charged hadrons down to ~30 MeV/c transverse momentum (pj) without 
the drawbacks from particles with small px curling up in the magnetic field. The nonzero px 
threshold is due to the 0.8 mm-thick beryllium beampipe, w^hich is not penetrable for pions and 
protons below px ~30 and 140 MeV/c, respectively. The loss of particles due to the beampipe 
is estimated to be less than 1 percent of the produced primary charged hadrons. 

Two analysis methods were used for the measurements of dN^h/dt] as a function of rj and 
centrality: one uses only pixel clusters in single BPIX layers (hit-counting method), and the 
other uses doublets of pixel clusters reconstructed from pairs of BPIX layers (tracklet method). 

The application of the pixel hit-counting method is a demonstration of the excellent pixel de- 
tector response and of its low occupancy even in this high-multiplicity environment, as well as 
the absence of noise and background. This method is not sensitive to detector misalignment 
or vertex-position resolution. The tracklet method is essentially a coincidence version of hit- 
counting. Using the angular coincidence of two hits from the same particle in different layers 
of the BPIX has the important feature of suppressing random noise. 



2 2 Trigger and event selection 

The paper is organized as follows. The triggering and event selection requirements are ex- 
plained in Section|2J followed by the description of the determination of the reaction centrality 
in Section|3J Sections|4]and|5]introduce the hit-counting and the tracklet methods, respectively. 
The systematic uncertainties are discussed in Section 16} while the final results are presented in 
Section [71 

2 Trigger and event selection 

The expected cross section for PbPb hadronic inelastic collisions at y^s^ = 2.76 TeV is 7.65 h, 
according to the chosen Glauber MC parameters described in Section |3J Electromagnetic inter- 
actions of the colliding nuclei at large impact parameter (ultraperipheral collisions, UPC) can 
lead to the breakup of one or both Pb nuclei with a much higher cross section. 

Minimum-bias (hadronic inelastic) collisions were selected by the Level-1 trigger system com- 
bining the logical OR of two clean and highly efficient triggers. One of them was the BSC 
coincidence, which requires at least one segment of the BSC firing on each side of the interac- 
tion point. The other was an HP coincidence trigger, which requires at least one HP tower on 
each side to have deposited energies that exceed the readout threshold. Both triggers accept 
noise at a low rate (less than 1 Hz with two noncolliding beams at full intensity), and have a 
very high efficiency (approximately 99% after the requirement of a reconstructed vertex). In 
order to suppress noncollision-related noise, cosmic-ray events, radioactivation, instrumental 
multiple triggering effects, and beam background, two colliding ion bunches were required to 
be present in coincidence with each one of these triggers, using information from the BPTX 
devices. The HP and BSC coincidence triggers were found to be largely insensitive to single- 
dissociation UPC, as discussed at the end of this section. 

The collision rate was 1.0-1.85 Hz per colliding bunch pair during the PbPb data taking period. 
Therefore, w^ith an orbit frequency of 11 245 Hz, the average number of collisions per bunch 
crossing was 0.9-1.6 x 10^^. There were 129 x 129 colliding bunches in the LHC at the time of 
data taking with no CMS magnetic field. 

In order to reject beam-gas interactions, large-hit-multiplicity beam background, and UPC, sev- 
eral offline event selection requirements w^ere imposed: 

• Events containing beam-halo muons and other particles from upstream collisions 
were identified and excluded from the analysis, by requiring the time difference 
between two hits from the BSC stations on opposite sides of the interaction point to 
be within 20 ns of the mean flight time between them (73 ns). 

• The large-multiplicity beam-background events were removed by requiring the com- 
patibility of the observed pixel-cluster lengths (defined in Section|4]| with the hypoth- 
esis of a PbPb interaction. This filter is the same as the one used in Ref. HI. 

• An HP coincidence requirement was imposed. At least 3 HP towers w^ere required 
on each side of the interaction point with at least 3 GeV total deposited energy in 
each tower. 

Purthermore, the presence of a reconstructed event vertex was required. The anal- 
ysis methods use their corresponding analysis objects (pixel clusters and tracklets, 
respectively) to reconstruct the interaction point. The vertex reconstruction is only 
done along the beamline; the transverse position of the vertex is taken to be that of 
the beam axis |J5l. The methods to determine the collision vertex are described in 
Sections E] and im 



• 



The measurement of the dNch/drj distributions was performed using 100 031 events, corre- 
sponding to an integrated luminosity of 13 mb^^. Correction factors were determined using 
simulated events generated with the ampt Monte Carlo (MC) [6J program. This program com- 
bines the HIJING event generator [7] with the ZPC parton cascade procedure [8| and the art 
relativistic transport model ||9l for the last stage of parton hadronization. The default tune is 
used, and the simulated events are reconstructed with the same version of software as used to 
process the collision data. This event generator produces a larger tail in the multiplicity distri- 
bution than that observed in data, making the entire observed multiplicity region completely 
covered in the simulation. The charged hadron multiplicity in the most-central collisions is 
20% higher in AMPT than in data, but since the analysis is done in bins of multiplicity, it is 
insensitive to this difference. 

The event selection for hadronic collisions was fully efficient for (mid)central PbPb collisions. 
For peripheral collisions the event selection efficiency was determined by comparing periph- 
eral PbPb data and a/s = 2.76 TeV pp data with the AMPT and PYTHIA Z2 QUI simulations. 
Based on these studies, the total event selection efficiency of the minimum-bias trigger for 
events produced in hadronic PbPb interactions was found to be (99 ± 1)%. 

The UPC contamination in the selected event sample was estimated using the photo-dissocia- 
tion simulations from Ref. [llj. The single-lead photo-dissociation events were found to be 
100% rejected by the event selection criteria outlined above, while half of the double-lead 
photo-dissociation events are found to pass the minimum-bias trigger. Such a UPC contam- 
ination amounts to (1 ± 0.5)% of the total number of events collected and populates ?»15% 
(5%) of the 95-100% (90-95%) most peripheral (largest-centrality) events, being negligible for 
the remaining 0-90% fraction of the PbPb cross section. 

3 Centrality determination 

In studies w^ith heavy ions, it is important to determine the degree of overlap of the two collid- 
ing nuclei, the so-called centrality of the interaction. Centrality is estimated using the sum of 
transverse energy in tow^ers from both HP at positive and negative z positions. The distribution 
of the total transverse energy, after the trigger efficiency and the UPC corrections, was used to 
divide the event sample into bins, each representing 5% of the total nucleus-nucleus interaction 
cross section. The bin corresponding to the most central events (i.e. smallest impact parameter) 
is the 0-5% bin, the next one is 5-10% and so on. The distribution of the HP signal, along with 
the cuts used to define the various event classes, is shown in Fig. IT] The UPC are concentrated 
in the two most-peripheral bins. To avoid them completely, only the 0-90% bins are used for 
the measurements reported in this paper. 

The centrality binning using equal fractions of the total interaction cross section can be corre- 
lated with more detailed properties of the collision. The quantity of interest for this measure- 
ment is the total number of nucleons in the tw^o Pb nuclei that experienced at least one inelastic 
collision, Npart- The average values of Npart for the various centrality bins (from most-central 
to most-peripheral), together with their uncertainties, are given in Table [l] The Npart values 
w^ere obtained using a Glauber MC simulation ||T2llT3l with the same parameters as in Ref. |[T4)|. 
These calculations were translated into reconstructed centrality bins using correlation func- 
tions between Npart and the measured total transverse energy, obtained from AMPT simulated 
events. Different Glauber MC samples were produced varying the Glauber parameters within 
the uncertainties from Refs. [15J and ||T6l . The variation in the final results is quoted as the 
uncertainty in Npart- 
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Figure 1: Distribution of the total transverse energy in the HF used to determine the centrality 
of the PbPb interactions. The centrality boundaries for each 5% centrality interval are shown 
by the dashed lines. 



Table 1: Average Npart values and their uncertainties for each PbPb centrality range defined 
in 5 percentile segments of the total inelastic cross section. The values were obtained using a 
Glauber MC simulation with the same parameters as in Ref. [I4]. 



Centrality 

Npart 


0-5% 
381 ±2 


5-10% 
329 ±3 


10-15% 
283 ±3 


15-20% 
240 ±3 


20-25% 
203 ±3 


25-30% 
171 ±3 


Centrality 

Npart 


30-35% 
142 ±3 


35-40% 
117 ±3 


40-45% 
95.8 ± 3.0 


45-50% 
76.8 ± 2.7 


50-55% 
60.4 ± 2.7 


55-60% 
46.7 ±2.3 


Centrality 

Npart 


60-65% 
35.3 ± 2.0 


65-70% 
25.8 ± 1.6 


70-75% 
18.5 ±1.2 


75-80% 
12.8 ± 0.9 


80-85% 
8.64 ± 0.56 


85-90% 
5.71 ± 0.24 



4 Hit-counting method and corrections 

Charged particles traversing the pixel detector deposit a certain energy in the silicon sensors, 
resulting in a proportional amount of charge collected in the pixel readout cells. Contiguous 
pixel cells with charge above the readout threshold are merged into a pixel cluster. A pixel clus- 
ter might be split into multiple clusters if one of its pixel cells fluctuates below the threshold. 
This phenomenon is called cluster splitting. The fraction of split clusters w^as estimated from 
the cluster-to-cluster distance distribution. The fractions in data and simulation were found to 
differ by less than 0.6%. The pixel-cluster efficiency (i.e. the probability that a cluster is detected 
once a charged particle crosses a pixel-detector layer), as well as the fraction of large clusters 
split into two are important quantities for the measurement. 

The pixel-cluster efficiency has been extensively studied in pp collisions IH [TtII , indicating an 
efficiency of (99.5 ± 0.5)%. Despite larger particle multiplicities in PbPb collisions, the occu- 
pancy of the pixel detector is still smaller than 1% owing to its high granularity (whereas in the 
strip detector it reaches 20%), and the pixel detector exhibits the same excellent perfomance in 
PbPb as in pp collisions. 

The hit-counting measurement method is based on the correlation between the cluster length 
in z and the pseudorapidity of a particle originating from the interaction point. It measures the 
primary charged hadron multiplicity distributions using the occupancy of a certain layer of the 
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Figure 2: Left: Distribution of the angle-corrected pixel-cluster charge in units of equivalent 
kilo-electrons from 2.76 TeV PbPb data and simulation. Right: Pixel-cluster length along the 
beam direction in units of pixel cells for hits from the first layer of the BPIX, as a function of tj 
after the event selection. The solid red line shows the selection on the minimum cluster length 
used in the analysis. 



pixel detector by counting the reconstructed hits. The hit-counting method gives three largely 
independent measurements for the three barrel layers. A similar method was used by the 
PHOBOS experiment at RHIC |18| and also by CMS for earlier pp analyses [4, 17]. One of the 
disadvantages of the method is the strong reliance on detector simulation for correction factors. 
Therefore, the detector simulation was extensively studied and carefully compared to data. The 
simulation was found to give a very good description of the data in all observables related to 
detector performance. Such an example can be seen in the left panel of Fig. [2l which shows 
the distribution of the pixel-cluster charge in the first BPIX layer from data and simulation. 
Clusters were selected according to the cluster-size selection described in Section |4.2[ and the 
cluster charge was normalised by the impact angle estimated from the cluster location and 
vertex position. The simulation describes the data well over six orders of magnitude. 

4.1 Primary vertex reconstruction using clusters 

There is a linear relationship between the length of a pixel-cluster along the beam direction 
and the z position of the cluster. Thus, hits from primary tracks leave a characteristic V-shaped 
pattern in the plane of cluster size versus z position. Nonprimary hits (e.g. due to secondary 
particles or nuclear interactions) fall mostly outside this V-shaped region. Thus, a V-shaped 
band is used to scan the z axis; the z position with the largest number of associated clusters 
is used as the vertex z position. The vertex z position is thus obtained by maximizing the 
consistency of the pixel-cluster lengths and global z positions with a primary vertex hypothesis. 

4.2 Cluster selection 

Particles travelling from the primary vertex at a small angle w^ith respect to the beam axis pro- 
duce larger clusters in the BPIX layers than those at large angles. The cluster size is proportional 
to I sinh;y|, where the pseudorapidity t] of the cluster is computed with respect to the recon- 
structed vertex (right panel of Fig. [2]|. Particles from background processes (decays in flight, 
nuclear interactions, etc.) often have smaller clusters than those produced in the primary in- 
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teraction, since their crossing angle is not correlated to the rj of the hit. Thus, a large fraction 
of clusters from background processes (refered to as background clusters) can be rejected by a 
selection based on the cluster size variable. The selection is defined in t] bins (shown as a red 
line in the right panel of Fig. |2]). 

4.3 Corrections 

Not all of the background clusters are removed by the cluster selection described above since 
they can occasionally mimic the length of clusters generated by primary particles. The correc- 
tion factor ;^(?/,M) is defined as the ratio of the number of selected clusters in the data to the 
number of primary charged hadrons at a pseudorapidity f] and a given cluster multiplicity M. 
It is calculated from simulation using: 

where M denotes the total number of clusters passing the cluster selection, N^f{ri,M) is 
the number of selected clusters, and ^hMy-oniV'^) the number of primary charged hadrons 
in the simulation. This correction factor is used to convert the measured Nhit(?/, M) pixel- 
cluster distributions from data into the corresponding primary charged hadron distributions, 

Nhadron(j/.M). 

The xiV' ^) correction factor is only weakly dependent on the physical process producing the 
hadrons, since it mainly contains information on the detector geometry. In a perfectly hermetic 
and 100% efficient detector, x{^, M) would be slightly above unity, as not only the primary, but 
also the secondary particles can generate hits. For detectors covering a limited solid angle, its 
values will be between and 1. For very large multiplicities, ;t(//, M) may decrease with in- 
creasing M because of the more significant occupancy (provided the primary/secondary ratio 
stays roughly constant). However, the occupancy of the silicon pixel layers is observed to be 
small and no apparent decrease oix{^, M) is visible with increasing centrality. The x correction 
increases with increasing distance from the interaction point, because layers further from the 
primary vertex are hit by more decay products, as well as by more secondaries from nuclear 
interactions. 

The correction factor in the first layer is in the range 1.0-1.2, while its average value in the 
0-10% centrality bin is 1.1, 1.2, and 1.3 for the first, second, and third layers, respectively. 

The pseudorapidity distribution of charged particles, for a fixed M, is calculated from the mea- 
sured Nhit('/, M) distribution, correcting for the hit/primary charged hadron ratio and normal- 
ising it to the number of events with multiplicity M passing the event selection, Nseiected(^)- 

^^(„ M) = ^ ^^^('^-^) (2) 

dri ^'' ' A^^(?/,M)N3eiected(M)' ^ ^ 

where A)/ is the width of the rj bin. 

The event selection efficiency for a given multiplicity is determined as the ratio of the number 
of MC events with multiplicity M which pass the event selection criteria ^^iected(^) ^^ ^^^ 
total number N^^{M) generated with multiplicity M: e(M) = Ng^iected(^)/^tot'^(^)- 

To measure the final, multiplicity-independent pseudorapidity distribution, the multiplicity- 
dependent distributions derived from Eq.^are weighted by the event selection efficiency e(M) 



and then summed over M: 



y^) ~ 5-- AT /I ,r\ 1 ■ ^'-'' 



dV EMNselected(M) 



e(M) 



Because a reconstructed event vertex is required as part of the event selection, the sum is over 
M>0. 

5 Tracklet method and corrections 

Tracklets are two-hit combinations in different layers of the BPIX that are consistent with a 
particle originating from the primary vertex. The tracklet analysis makes use of the correlation 
between hit positions: pairs of hits produced by the same charged particle have only small 
differences in the pseudorapidity (A?/) and the azimuthal angle (A(^) with respect to the primary 
vertex. 

5.1 Primary vertex reconstruction using tracldets 

In this method a tracklet-based vertex finder is used. In the first step, a hit from the first BPIX 
layer is selected and a matching hit is sought. If the magnitude of the difference in azimuthal 
angle (Ai^) between the two hits is smaller than 0.08, the pair is saved as a proto-tracklet. This 
procedure is repeated for each first-layer hit to get a collection of proto-tracklets. For each 
proto-tracklet, the expected longitudinal vertex position is found using: 

z = zi -ri(z2-zi)/(r2-ri), (4) 

where Zi(2) is the z position of the first (second) layer hit, and ri(2) is its radius. The calculated 
z positions are saved as vertex candidates. The second step is to determine the primary vertex 
from the vertex candidates. If the magnitude of the difference between the z positions of any 
two vertex candidates is less than 0.14 cm, they are combined as a vertex candidate cluster. 
Finally, the vertex candidate cluster with the highest number of vertex candidates is selected as 
the primary vertex. The final vertex z position is determined by the average z position of the 
vertex candidates in the cluster. 

5.2 Tracl<let reconstruction 

All three barrel layers of the pixel detector are used in pairs: lst+2nd, lst+3rd, and 2nd+3rd. 
The differences in pseudorapidity and azimuthal angle, as w^ell as the two-dimensional separa- 
tion AR = ^J[l^r]Y + {^<PY betw^een the two hits of a tracklet, are important for characterising 
the tracklet. 

Tracklets are reconstructed in three steps: 

1. For each reconstructed hit, the pseudorapidity is calculated using the primary vertex lo- 



cation. Hits that pass the cluster size selection (as described in Section 4.2 1 are kept for 
further analysis. 

2. Starting with a reconstructed hit in the fl* layer and looping over the reconstructed hits 
in the &* layer (with h > a), all possible combinations with |AR| < 0.5 are saved as 
proto-tracklets. 
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3. Proto-tracklets are sorted in AR. If a &*-layer hit is matched more than once, the proto- 
tracklet with the smallest AR is kept. The selected proto-tracklets are the final recon- 
structed tracklets. 

In addition to primary charged particles, the set of tracklets also include contributions from 
secondary interactions in the beampipe, particles from weak decays, and combinatorial back- 
ground. 

The combinatorial background tracklets are defined as combinations from secondary hits and 
hits from different primary tracks. The background fraction is largely suppressed by the AR 
ordering and the selection of tracklets (described in the next section). The tracklets from sec- 
ondary particles are suppressed, and the correction for the remaining contribution relies on 
simulation. 

5.3 Combinatorial and secondary particle background 

Background tracklets can be created from incorrectly associated hits. The Arj and A^ of a track- 
let are very useful quantities for the separation of signal and combinatorial background track- 
lets. Because of the absence of magnetic field, selecting the best proto-tracklet with the smallest 
AR provides a powerful way to reject combinatorial background: signal proto-tracklets exhibit 
a correlation peak around AR = 0, while the background component extends to large AR. 

The Afj and A(p distributions of selected tracklets from minimum-bias collisions in data and 
simulation are shown in Fig.lSlfor combinations in the first and second pixel layers. The signal 
peaks at Arj and A^ = are clearly visible. Data and simulation show^ agreement over several 
orders of magnitude. 
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Figure 3: The (left) Arj and {right) Acp distributions for reconstructed tracklets in minimum-bias 
collisions from the first and the second pixel layers in data and simulation. 



The effect of secondary hits on the agreement between data and simulation seen in Fig. [3] 
was tested by adding random hits to simulated events. The simulated tracklet spectra w^ere 
found to be distorted even by a few^ percent of random hits, spoiling the agreement between 
data and simulation. Given the very good agreement in the tracklet spectra, the fraction /5 of 
combinatorial background tracklets can therefore be reliably obtained from simulation. The 
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value of fi is in the range 0-15%, depending on the multiplicity of the event, the pseudora- 
pidity of the tracklets, and the z position of the event vertex. The number of background- 
subtracted tracklets Mrackiet is determined from the raw number of tracklets in data NJj.^™j,jgj 
using: Nt,ackiet = (1 - /5) X N-^^1^^. 

5.4 Efficiency and acceptance correction 

To calculate the number of hadrons from the number of tracklets, an efficiency correction must 
be applied. The correction factor a(M, ?/,Zv) for the tracklet reconstruction efficiency is defined 
as 

^ ''"^ " [l-/3(M,,,z.)]N-:i^,f(M,,,z.)' ^'^ 

where Zv is the z position of the vertex, N^^^^^j^ (M, rj, Zy ) is the true number of charged hadrons 
in the simulated sample, and N^^^^^^^, is the raw number of selected tracklets in the MC sam- 
ple. The factor ix.{M,f],Zy) is used to calculate the charged hadron spectra from the measured 
background-subtracted tracklets. Typical values of a are less than 1.15 because of the high hit- 
reconstruction efficiency in the BPIX. At larger pseudorapidity, the correction factor increases, 
owning to the reduced acceptance. The size of the acceptance correction also depends on the 
position of the primary vertex. 

The pseudorapidity distribution of charged hadrons for a given multiplicity M is obtained 
from the measured number of tracklets (Nj^^ciet)' *^^ background fraction (/5), the efficiency 
and acceptance correction (a.), and the normalisation to the number of selected events: 



df] ^^' ^ A^ N3eiected(M) ' ^^^ 

where A^ is the width of the tj bin and Nsgiected(^) is the number of selected events used in 
each multiplicity bin. The a(l — /S) correction has a typical value larger than 0.85. For the final 
dNch/dr] distribution the multiplicity-dependent results are weighted by the event selection 
efficiency e{M) and then summed as in Eq.lsj 

6 Systematic uncertainties 

A summary of systematic uncertainties affected the measurement of dN^h/df] for the two anal- 
ysis methods is given in Table [2| 

The results from different BPIX layers (or layer combinations) from either of the measurement 
methods differ by less than 2%, and thus they are averaged using the arithmetic mean. The 
uncertainties of the averaged results from the hit-counting and tracklet methods are dominately 
systematic and largely correlated. Therefore, the two results are averaged using equal weights. 
Since the difference betw^een the measurements from the tw^o methods is smaller than 1%, the 
w^eighting procedure has very little effect on the numerical value of the final result. 

The uncertainties on the final average result are computed as follows. All the systematic uncer- 
tainties listed in Table |2] are correlated between the two methods, except those associated with 
the efficiency of reconstruction and misalignment. The correlated, (sh.c.); and (strackiet)// arid 
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the uncorrelated, (c^h.c.)/ srid (tT'trackiet); uncertainties of the hit-counting and tracklet methods 
(indexed by ;') are summed in quadrature: 



S = Y'S;[(Sh.c.); + (Stracklet);]^/2 and CT = y' (S;[(t7h.c.)f + (c^tracklet)f ]/2, (7) 

resulting in s and a correlated and uncorrelated uncertainties of the average, respectively. The 
total systematic uncertainty is then o'tot = V^^ + s^- 

In this paper the distributions of three observables are reported: dNch/dri\,j=o as a function 
of centrality class, (cfNch/d'/)/(Npart/2) as a function of //, and {dNch/ drjl^^o) / {Npari/2) as a 
function of Npart- 

The systematic uncertainties affecting the slope and those affecting the absolute scale of 
dNch/df]\,j=o and ((iNch/d?/|);=o)/(Npart/2) measurements are determined separately. System- 
atic uncertainty sources affecting the slope are those on the centrality and the Glauber calcula- 
tion of Npart; all other sources affect the scale. The results are presented in Section [7| with these 
two uncertainties shown separately. 

The slope of dN^hf dv\']=Q as a function of centrality is only affected by the uncertainty on 
the determination of the centrality bins. Both the slope and the absolute scale of the Npart- 
normalised distributions are affected by the uncertainty on Npart from the Glauber calculation, 
given in Table [Tj These contributions are computed by transforming the uncertainty in Npart 
into an uncertainty on the Npart-riormalised hadron density distribution using the derivative of 
the measured (dNch/rf?/)/ (Npart/2) distributions as a function of Npart- 

Table 2: Summary of systematic uncertainties on the dNch/dt] measurements and their sum for 
the two analysis methods. 



Source 


Hit-counting [%] 


Tracklet [%] 


Centrality (0-5% to 85-90%) 


0.5-15.6 


0.5-15.6 


Pixel hit efficiency 


0.5 


1.0 


Tracklet and cluster selection 


3.0 


0.5 


Acceptance uncertainty 


1.5 


1.5 


Correction for secondary particles 


2.0 


1.0 


Pixel-cluster splitting 


1.0 


0.4 


Reconstruction efficiency 


- 


1.9 


Misalignment 


- 


1.0 


Random hits 


1.0 


0.2 


Total uncorrelated uncertainties 


- 


2.1 


Total uncertainties 


4.2-16.2 


3.1-15.9 



The systematic uncertainties affecting the measurements are as follows. 

• Centrality: The determination of the centrality bins as a percentage of the total hadronic 
cross section relies on the hadronic event selection efficiency (99 ± 1)%, UPC contamination 
(1 ± 0.5)%, and the percentile binning of the centrality variable. Thus, the uncertainty in the 
event selection cross section causes uncertainty in the centrality binning by moving the bin 
boundaries, shifting the event population in each centrality bin. The effect of this centrality 
uncertainty on the final results was studied by repeating the analysis using various centrality 
tables (derived from the various trigger efficiencies allowed by the uncertainty on the effi- 
ciency). These studies indicate that the uncertainty of the dNch/df] values ranges from 0.5% for 
the 0-5% centrality bin to 15.6% for the 85-90% centrality bin. 
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• Pixel hit efficiency: The efficiency of the pixel layers is larger than 99%, measured from pp 
data II17|| . The pixel detector has low occupancy even in central heavy-ion collisions because 
of its fine segmentation. Therefore, the same pixel hit efficiency and uncertainty measured in 
pp collisions are used here. The pixel hit efficiency affects tracklets more, since two layers are 
required; a 0.5% inefficiency per pixel layer leads to a 1% inefficiency for tracklet finding. 

• Tracklet and cluster selection: Varying the cluster selection requirements (pixel-cluster length 
selection) and tracklets selection (requirement on AR) is used to estimate the uncertainty due 
to cluster and tracklet selection. The cluster selections were changed by one pixel unit in all tj 
bins and the AR selection by a factor of three. The observed differences in the final results (3% 
and 0.5%, respectively) are quoted as conservative systematic uncertainties. 

• Acceptance uncertainty: The positions of the BPIX modules are only slightly different in 
data and in simulation, but hits at the extreme edges of the BPIX are not used in the analysis, 
limiting the systematic uncertainty from this effect. The rj, (p acceptance was estimated from 
data by using the endpoints of tracklets to map the active surface of the BPIX layers. From this 
study, the acceptance correction is estimated to be 1% in pp collisions. In PbPb collisions (due 
to the longer luminous region in the beam direction) this uncertainty was increased to 1.5%. 
No correction is applied, but the effect is included in the systematic uncertainty. 

• Corrections due to hits from secondary particles: The sensitivity of the correction factors 
applied to remove hits caused by secondary particles was tested using two largely different 
event generators: AMPT and hydjet |19| . The relative fraction of strange particle production 
differs by 60% in the two generators, but the effect on the correction factor was found to be only 
2% for the case of the hit-counting analysis. The tracklet analysis is less sensitive to secondaries. 

• Pixel cluster splitting: The relative fraction of split clusters was estimated from the cluster- 
cluster distance distribution. This study shows that the number of split clusters in data does 
not exceed that in simulation by more than 0.5-0.7%. No correction is applied for this effect 
in the analyses, but a conservative systematic uncertainty (1% and 0.4%, respectively, for the 
hit-counting and tracklet analyses) is assigned. 

• Efficiency of tracklet reconstruction: The uncertainties in the simulation of the px and mul- 
tiplicity (M) distributions influence the determination of the tracklet reconstruction efficiency. 
The uncertainty (1.9%) is estimated based on variations of these quantities within reasonable 
limits: (pj) 'was modified by 10%, and the multiplicity variable was changed from using clus- 
ters to HF towers. 

• Misalignment: The hit-counting method is not sensitive to detector misalignments. The 
tracklet method has a sensitivity through the AR selection, which was studied by moving the 
reconstructed hit positions (the entire detector) by 0.3 mm in the simulation, while keeping the 
vertex position at the same place, giving a 1% change in the final result. Since the Ai^ and Atj 
distributions and the correlation widths agree well, no significant misalignment is seen in the 
data. 

• Random hits: With the restrictive event selection criteria, the contamination from beam-gas 
(high-occupancy) events in the final data sample is negligible. The other potential source of 
background is the accidental overlap between beam-gas and PbPb collisions. A conservative 
systematic uncertainty of 1% is assigned for the hit-counting analysis, which is more sensitive 
to this overlap than the tracklet analysis, for which 0.2% is assigned. 
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Figure 4: Le/it: dNch/dri\,j=o as a function of centrality class in 2.76 TeV PbPb collisions from this 
experiment (solid circles) and from ALICE (open squares) |20|. The inner green band shows the 
measurement uncertainties affecting the scale of the measured distribution from this analysis, 
while the outer grey band shows the full systematic uncertainty, i.e. affecting both the scale and 
the slope. Right: Measured dNeh/rf?// (Npart/2) distributions from this analysis as a function of 
rj in various centrality bins. 



7 Results 

The hit-counting and tracklet dN^.-^/ dtj results are in good agreement; their average difference 
is smaller than 1%. Their individual results are averaged as described in Section |6J and these 
averages are presented as the final results. 



The left panel of Fig. El presents the measured dN^i^/dr]\^=Q values as a function of centrality. 
The statistical uncertainties are negligible, while the systematic uncertainties are shown as two 
bands. The inner green band show^s the measurement uncertainties affecting the scale of the 
measured distribution, while the outer grey band shows the full systematic uncertainty, i.e. 
affecting both the scale and the slope. Details on the calculation of the uncertainty bands are 
given in Section[6] The charged hadron density for the 5% most-central events (0-5% centrality 
bin) is measured to be dNch/^?/ 1)^=0 = 1612 ± 55 {syst.). These results are consistent with 
those of ALICE [20J within the uncertainties, as shown in Fig. [4] (left). The error bars of the 
ALICE points in the figure show the total statistical and systematic uncertainties. The CMS 
measurements cover the centrality range of 0-90%, extending the ALICE results (0-80%) to 
more-peripheral collisions. 

In order to compare bulk particle production for different colliding nuclei and at different en- 
ergies, the charged-hadron density is divided by the average number of participating nucleon 
pairs, Npart/2, determined for each centrality bin. The Npart values are obtained using the 
Glauber calculation, by classifying events according to their impact parameter, without refer- 
ence to a specific particle production model (Table [l]). 

The measured {dN^h/ dt]) / {Npart/'^) distributions as a function of rj in various centrality bins 
are shown in the right panel of Fig.|4] The uncertainty bands of these distributions also include 
the Glauber uncertainty on Npart- The )] dependence of the results is weak, varying by less than 
10% over the t] range. The slight dip at t] = is a trivial kinematic effect (Jacobian) owing to 
the use of pseudorapidity (rj) rather than rapidity (y). 
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Figure 5: Left: Measured {dNch/drj\^^o) / {Npart/'^-) as a function of the number of participants 
in 2.76 TeV PbPb collisions from this analysis and the ALICE experiment [20], from RHIC 
im at 200 GeV and 19.6 GeV, and from extrapolated pp results from CMS |17| and ALICE 
1122 1 . Systematic uncertainties affecting the scale of the measurements from this analysis are 
shown as inner green error bands and the total systematic uncertainties as an outer grey band, 
while the error bars indicate statistical uncertainties. The black stars are shifted slightly to 
the right for better visibility. The ALICE and the averaged RHIC results are from [20| and 
II2TI , respectively. Right: Results from this analysis are compared with model predictions of 
{dNch/ df]\,j=o) / {Npait/'^) as a function of the number of participants in 2.76 TeV PbPb colli- 
sions. The model predictions are taken from Refs. p3| , 124), and 



The left panel of Fig. [S] presents the measured {dN^h/ dr]\ 1^=0) / {^pan/ '2-) as a function of Npart- 
The statistical uncertainties on the CMS results are indicated by error bars (negligible), while 
the systematic uncertainties are shown as two bands. The inner green band shows the system- 
atic uncertainty affecting the scale, while the outer grey band shows the full systematic uncer- 
tainty. The error bars on the ALICE [201 and the RHIC ||2T) points show the quadratic sum of 
the statistical and systematic uncertainties. The RHIC results are multiplied by numerical fac- 
tors to match the Npart-normalised multiplicity observed at the LHC for central collisions. The 
pp results shown in the figure are obtained from the measured non-single-diffractive (NSD) 
dNa,/df]\^=o = 4.47 ±0.2 (CMS) [17] and the inelastic dNa,/df]\^=o = 3.77+°;^^ (ALICE) EH 
values at 2.36 TeV, using the ^ys dependence of the charged hadron multiplicity density mea- 
sured in NSD and inelastic collisions from Ref. ID. The error bars on the pp points show the 
total (statistical and systematic) uncertainties. The Npart values used for the normalisation by 
CMS and ALICE differ by less than 2%. Within the uncertainties, the Npart-normalised hadron 
densities follow a similar dependence on centrality for all centre-of-mass energies, although 
the lower-energy collider data appear to have a flatter dependence on Npart- 

The phenomenological descriptions of particle production in nuclear collisions are often based 
on two-component models, combining contributions from perturbative QCD processes, i.e. 
(mini)jet fragmentation and soft interactions. The data are compared to three different ap- 
proaches: (i) HIJING 2.0 [23], which basically scales (via the number of incoherent nucleon- 
nucleon collisions) the (semi)hard parton scatterings and fragmentation (Lund model [26]) im- 
plemented in PYTHIA after accounting for the "shadow^ing" of the nuclear parton distribution 
functions; (ii) parton saturation approaches ||24t , which model heavy-ion interactions as the 
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collision of two dense multigluon wavefunctions with cross sections peaking at a semihard 
scale (saturation momentum of ?»2-3 GeV/c at the LHC) I|27ll28| , followed by their fragmenta- 
tion according to a simple parton-to-hadron local-duality prescription; and (iii) the DPMJET-III 
MC program 1125 1 , based on the Regge-Gribov theory. This is an extension of the phojet 1129 1 
program in w^hich interactions from soft degrees of freedom (Pomerons) can fuse in the dense 
initial state. They are extended consistently into the hard regime via "hard" or "cut" Pomerons, 
and then fragmented using the standard Lund model. 

The measured {dNch/ drj\,j^Q) / {Npait/'^) versus Npart distribution is compared to the various 
model predictions in the right panel of Fig. [5] The two-component HIJING 2.0 model, w^hich 
has been tuned to high-energy pp and central PbPb data, describes the general shape of the 
data. The HijiNG model includes an impact-parameter-dependent gluon shadowing parame- 
ter s^, which limits the rise of particle production with centrality The magnitude of the particle 
production favours a relatively large Sg = 0.23 value, although the shape of the centrality 
dependence prefers a smaller Sg = 0.17. The observed centrality dependence is well repro- 
duced by the saturation model of Ref . [24J . Both Refs. 1123 1 and p4) were published knowing 
the result of ALICE [30] on the multiplicity of the 5% most-central collisions, although pre- 
vious saturation-based calculations (e.g. [31J ) predicted central charged hadron densities very 
similar to those finally measured. The DPMJET-III model is capable of describing the charged 
hadron multiplicity in the most-central collisions, but shows a stronger rise with centrality than 
observed in the data. The measured particle densities provide basic constraints on the initial 
conditions of the quark-gluon plasma in any hydrodynamical approach employed to study 
PbPb interactions at the LHC 
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Figure 6: Normalized charged hadron pseudorapidity density at ?/ = as a function of centre- 
of-mass energy for the 0-5% most-central nucleus-nucleus collisions (top set of points) and 
0-70% centrality (middle set), and for NSD pp collisions (bottom set). The fits to power-law 
functions are shown by lines. The grey band around the PbPb CMS points indicates the total 
systematic uncertainty. The statistical uncertainty is negligible. The error bars on the other 
points indicate statistical and systematic errors. The ALICE, PHENIX, and PHOBOS results 
(w^hich are shifted slightly to the right for better visibility) are taken from Refs. [l20l , l|2T1 . and 
1133 J , respectively The NSD pp results of CMS, ALICE, UA5, and UAl are from Refs. EUlZl, 
1, 134.1 , and Il55l , respectively. 



The collision-energy dependence of the measured (dNch/t^?/|);=o)/(Npart/2) for 0-5% and 0- 
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70% centrality from this analysis and from ALICE, PHENIX, and PHOBOS can be seen in 
Fig. |6] The PHENIX and PHOBOS points are taken from Refs. [21J and ESI, their error bars 
representing both the statistical and systematic uncertainties. Systematic uncertainties of the 
measurements from this analysis are shown as an error band, while the statistical uncertain- 
ties are negligible. The NSD pp results of CMS, ALICE, UA5, and UAl are from Refs. BIITtI. 
f22\ , 1341 , and [35 1, respectively. The Npart values at different collision energies are different 
for a fixed centrality bin. When the Npart dependence of {dNch/ drj) / {Npart/2) from PHENIX 
and PHOBOS are used to extrapolate their (dNch/'i?/)/(Npart/2) results shown in Fig.|6]to the 
Npart values appropriate for the LHC, they change by no more than 3%. This correction is not 
applied in Fig. |6] The normalised charged hadron densities shown in Fig. [6] are fit to a power- 
law function: a + s" . The fit returns the value n = 0.13 for PbPb and n = 0.10 for NSD pp 
collisions. The results of the fits are shown by the straight lines in Fig. pj These results provide 
additional constraints on the energy evolution of the saturation momentum in the proton and 
nuclei ||27l [28| , as well as in general on the pj cutoff between soft and hard dynamics used in 
the models of particle production in high-energy hadronic collisions. 

8 Summary 

A measurement of charged hadron multiplicity as a function of pseudorapidity and centrality 
in PbPb collisions at ^/s^ = 2.76 TeV has been reported. For the 5% most-central collisions, 
a primary charged hadron density of 1612 ± 55 is measured, which represents an increase of 
a factor of 3 compared to similar measurements at RHIC energies. The dNch/dt] distributions, 
measured over the range |?/| < 2.5, show^ w^eak rj dependence, the variation being less than 
10%. The Npart-normalised multiplicity distributions from RHIC and the LHC have a similar 
dependence on centrality, although the lower-energy collider data has a somewhat flatter de- 
pendence. A parton saturation model describes well the observed centrality dependence. The 
collision-energy dependence of the measured hadron multiplicities at central rapidities is well 
modelled by a power-law function of the type a + s'^^ . These results provide information on 
the parton structure of the nucleus and the proton and its evolution as a function of centre-of- 
mass energy. They also give additional constraints on the initial conditions in nucleus-nucleus 
collisions at LHC energies for hydrodynamical evolution studies of the strongly interacting 
produced system. 
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